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We analyse deformation mechanism in neutron-rich Cr, Fe and Ti isotopes with iV=32- 
44 by means of a Skyrme-Hartree-Fock-Bogoliubov mean-field code employing a two- 
dimensional mesh representation in the cylindrical coordinate system. Evaluating system- 
atically the quadrupole deformation energy, we show that the Skyrme parameter set SkM* 
gives a quadrupole instability around the neutron number iV~38-42 in Cr isotopes, where 
the deformation energy curve suggests a transitional behavior with a shallow minimum ex- 
tending to a large prolate deformation. Roles of a deformed -/V=38 gap and the position of 
the neutron go orbit are analysed in detail. 

P5 , §1. Introduction 



Deformation in the nuclear shape is one of the fundamental features of the 
nuclear structure. Mechanism of the deformation reflects quantal effects such as the 
shell effect, the pairing correlation, and the long range proton-neutron correlation. 
. Thus, analyses of a new deformation region often provide us with new insights into 

nuclear structure. A recent example is the quadrupole deformation in neutron-rich 
Z^12 isotopes, which have been a crucial probe to disclose the weakening of the 
iV=20 shell gap in the neutron-rich nuclei. 1 ) More recently another possibility of a 
\ new deformation region has been suggested in neutron-rich Cr isotopes 2 ) - . 7 ) 

The self-consistent mean-field theory is a powerful scheme to analyse the defor- 
mation mechanism when we deal with nuclei far from the stability line in medium and 
heavy mass regions. 8 ) There are three important features which should be taken into 
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account in this case: the pairing correlation, the continuum states and the deforma- 
tion. The pairing correlation can be described by using the Bogoliubov's generalized 
quasi-particle scheme while the conventional BCS approximation is not suitable for 
nuclei far from the stability line. The continuum states become important in nuclei 
near the drip-lines where the Fermi energy becomes close to zero. In this case one has 
to treat continuum quasi-particle states which have spatially extended wave func- 
tions. These two requirements are fulfilled in the Hartree-Fock-Bogoliubov(HFB) 
formalism once it is represented in the coordinate-space. 9 )> 10 ) To describe the defor- 
mation, then, the HFB method has to be extended to a two- or three-dimensional 
problem, for which several techniques using the truncated Hartree-Fock basis, n )~ 13 ) 
the transformed harmonic oscillator basis, 14 )> 15 ) and the B-Spline basis 16 )' 17 ) have 
been developed recently. 
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In the present paper, we develop our own 2D coordinate-space Skyrme-HFB 
code, which enables us to describe axially deformed unstable nuclei. The formal- 
ism is similar to that of Ref. 16), which utilizes the B-Spline and Galerkin methods 
on top of the cylindrical coordinate system. In the present approach, we adopt a 
two-dimensional mesh representation using the cylindrical coordinate system and a 
multi-points formula for differential operators. Since the 2D mesh approach using 
the cylindrical coordinate has not been explored except a few Hartree-Fock 18 ) and 
Hartree-Fock-Bogoliubov calculations, 19 )> 20 ) we shall test numerical accuracy of our 
code. Then, as an application of this code, we shall investigate the possible deforma- 
tion region around the neutron-rich Cr isotopes. A global calculation of the ground 
state deformation including the neutron-rich Cr isotopes is performed by means of 
the Skyrme-HFB 15 ^' 21 )' 22 ) and the relativistic mean-field models. 23 ) Here we perform 
detailed analysis of the deformation mechanism by looking into the deformation en- 
ergy curve and the shell structure of the single-particle orbits, and dependence on 
the Skyrme interactions. 



2.1. The HFB equation in the cylindrical coordinates system 

The HFB equation in the coordinate space representation is written as 



neutrons and protons (q = n or p) with a and a' being spin indexes (a, a' = |,|). 
In the present paper we assume the axial symmetry of nuclear deformation and 
the z-axis is chosen as the symmetry axis. The quasi-particle wave function has a 
quantum number Q for the z-component of the total angular momentum. Using the 
cylindrical coordinate (r, z, tp), the quasi-particle wave function of the n-th eigenstate 
is written as 



§2. The HFB formalism 





(2-2) 



Expressing the Hearted-Fock 

h(r, z, if) as 

h g (r,z,tp) ■- 



Hamiltonian h(r, z, ip) and the 



pairing Hamiltonian 



e^h\Jr,z) 



e- il ?h q n (r,z) 



(2-3) 



h q (r,z,ip) 



e^h\Jr,z) 



e-^h q n (r,z) 
~ h l\( r > z ) 



(2-4) 
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the HFB equation in the cylindrical coordinate representation is written as 

H q ^ n Q q = E nQq $ nQq (2-5) 



where 




We adopt the Skyrme force and h q , (r, z) is expressed 16 ' in terms of the force param- 
eters and the normal density p njP (r,z), the kinetic density T n>p (r, z), the spin-orbit 
density J n ,p(f, z). We adopt the density dependent delta interaction(DDDI) 



v P air = \v (l - P.) (l - (^V) 8(r - r') 

for the effective interaction producing the pairing potential h q , , where h 9 acr , is ex- 
pressed in terms of the total nucleon density p(r,z) and the pair density p q (r,z). 
The eigenstates and eigenenergies of the quasi-particles are obtained by solving the 
HFB equation (2-5) for each 17. Since we assume the time reversal symmetry, we 
need to solve only for 1? > 0. 

We consider a rectangular area < r < r max and — z max < z < z max in the 
(r,z) plane. The r-axis is discretized at r = h/2, 3h/2, ■ ■ ■ ,r max with N r mesh 
points, while the z-£ixis is discretized from — Zmax — -N z h to 

z zmax — N z h with the 

same equidistant interval h. The quasi-particle wave functions (f)(r, z)(the indices 
are omitted here for simplicity) and the HFB Hamiltonian H q are represented on 
this two-dimensional mesh. At the outer boundary z = ±z max and r = r max , 
we impose (f)(r max ,z) = and <p(r, ±z max ) = 0. At the origin r = 0, we impose 
4>(—r, z) = (—) lz (f)(r, z) determined by the parity (— ) lz where l z is the z-component of 
the orbital angular momentum. Note that we utilize shifted mesh points r = (n+^)h 
so that we do not need to evaluate the value of (f> at r = 0. We made a comparison 
with calculations where the r-axis is discretized on a slightly different mesh points r = 
: r max and we found that the adopted discretization with r — (n-\-^)h has 
better accuracy. We adopt the 11-points formula to represent differential operators. 
Concerning the integrals, we adopt trapezoidal formula since its is known that the 
high order formula does not always guarantee a better accuracy. 24 -* Consequently a 
quasi-particle wave function is expressed as a vector with N = 4N z N r dimension and 
the HFB Hamiltonian is a matrix having the size of iV 2 . We obtain the eigenvalues 
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and the eigenfunctions by diagonalizing this Hamiltonian matrix. The Hamiltonian is 
non-symmetric because of the finite-points representation of the differential operators 
around r = 0. We diagonalize the Hamiltonian using the QR method to solve the 
non-symmetric eigenvalue problem. 

Since we use the delta type pairing force, we introduce a cut-off with respect to 
the quasi-particle energy. The adopted value of cut-off energy is 60 MeV. We need 
also a cut-off with respect to the azimuthal quantum number 17, and we adopted 
the cut-off value of J? = In order to obtain the self-consistent solution, we utilize 
an iteration scheme. As an initial condition, we start with single-particle orbits in 
a deformed Woods-Saxon potential. As an initial pairing potential, we chose one 
having a Woods-Saxon shape. We revise only a fraction of the densities in each step 
of iteration, and this fraction is randomized within an interval 0.2-0.4. When about 
200 iterations pass, the increment of the energy reaches 1 x 10 -5 MeV, at which 
we stop the iteration. In order to obtain a deformation energy, we perform also a 
constrained HFB calculation. As a constraint operator, we use the mass quadrupole 
operator Q2, which is defined as 

Q 2 = y^ 2 y 20 (f) = 2z 2 -r 2 

{R = ^r 2 + z 2 ). We use the quadratic constraint method. > 

2.2. Test of the numerical code 

To check the accuracy of the code, we compare our result with those in the two 
previous Skyrme-HFB calculations. One is the THO basis method developed by 
Stoitsov et al. 14 ) and the other is the B-Spline basis method developed by Teran et 
al. 16 ) We choose neutron-rich Zr isotopes. Along with Refs.14) and 16), we adopt 
the Skyrme parameter SLy4 and the volume type pairing force, i.e. rj = 0. The box 
size is 12 fm for both r max and z max , and the mesh interval h is 0.6 fm. The pairing 
strength Vq is the same in Refs. 16) and 17). The obtained ground state properties 
of neutron-rich Zr isotopes are shown in Table I and Fig.l. The density and the pair 

Table I. The total binding energy BE, the Fermi energy X q , the average pairing gap A q and the 
RMS radius \J (r 2 ) in 102 Zr obtained in the present work. The values 16 ' of the B-Spline and the 
THO methods are compared. 





present work 


B-Spline 


THO 


BE [MeV] 


-856.92 


-859.61 


-859.40 


A„ [MeV] 


-5.42 


-5.46 


-5.42 


A p [MeV] 


-12.07 


-12.08 


-12.10 


A n [MeV] 


0.23 


0.31 


0.56 


A p [MeV] 


0.35 


0.34 


0.62 


vV> M 


4.58 


4.58 


4.58 



density for neutrons are shown in Fig.l. It is seen that the neutron density has a 
large quadrupole deformation and the neutron pair density is also deformed. Note 
that the pair density exhibits a complex structure in the inner region. 
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Fig. 1. (a) The density p n (r, z) and (b) the pair density p n (r, z) for neutrons in 102 Zr. The unit in 
the legend is fm~ 3 . 



In Table I, we compare our result for 1U2 Zr with the results of the B-Spline and 
the THO methods in Ref. 16). It is seen that our code reproduces reasonably well. 
Especially, the root mean square (RMS) radius \J (r 2 ) and the Fermi energies X njP 
are in a good agreement. Note however that there is a deficiency in the total binding 
energy by about 2-3 MeV. This is probably due to use of the moderate mesh size 
(h=0.6 fm) and the finite-points formula for the differential operators as a deficiency 
due to the finite mesh size has been known in previous HF+BCS calculations. 26 )' 27 ) 
There is some deviation in the average neutron pairing gap A n . It should be noted 
that the pairing gap is rather small ~ 300 keV. Generally the value of the pairing 
gap is sensitive to minor differences, especially when the pairing gap is small. A 
small difference in defining the cut-off energy may be suspected as we adopt the 
quasi-particle energy cut-off while the THO and B-Spline methods use so called the 
equivalent single particle energy cut-off. The finite mesh size effect could be another 
origin. 

To make a systematic comparison, we have calculated an isotopic chain of zir- 
conium for iV=102-122. The mass quadrupole moment Q2, the RMS radius, the 
neutron and proton pairing gaps, and the two-neutron separation energy are com- 
pared with those of Ref. 17) in Fig.2(a)-(d). In this calculation, the pairing strength 
is Vo=-187.1305 MeV fm~ 3 . The quadrupole moment (Fig.2(a)) as well as the RMS 
radius (Fig. 2(b)) agree very well with those in the B-Spline and the THO methods 
in the whole region (^4=102-122) shown in the figure. Especially the deformation 
is described well. The difference in the neutron separation energy is smaller than 
the 2-3 MeV difference in the absolute value of the total binding energy. The better 
accuracy in the relative quantities is consistent with the same feature found in the 
previous Skyrme-HF + BCS calculations using the 3D Cartesian mesh representa- 
tion. 28 ) The average pairing gap is also reproduced reasonably well as shown in 
Fig.2(c) although a small but non-negligible difference is noted. 

In Figure 3, we display the deformation energy curve as a function of the 
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A A 

Fig. 2. The comparison of (a) the mass quadrupole moment Q2, (b) the RMS radius a/ (r 2 ), (c) 
the neutron average pairing gap A n , and (d) the two neutron separation energy S^n in the 
neutron-rich Zr isotopes between the present work and the B-Spline and the THO methods. 17 ^ 

quadrupole deformation parameter (3 defined by (3 = \/f" )A • dependence 
of the quadrupole deformation energy on the mesh size is seen from the comparison 
of the results obtained with the mesh sizes 0.6 and 0.8 fm. If we shift the total energy 
by about 1.8 MeV, we see that the difference between the two curves becomes within 
about 200 keV. We thus conclude that the deformation energy can be evaluated to 
the accuracy of this order. 




-527 1 1 1 1 1 1 1 1 1 <— 

-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 
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Fig. 3. The deformation energy curves as a function of the deformation parameter (5 calculated for 
62 Cr using SkM*. The solid line is the case of h =0.6 fm and the dash line for 0.8 fm. 
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§3. Quadrupole deformation around the neutron-rich Cr isotopes with 

N ~38 



3.1. The neutron-rich Cr isotopes and neutron single-particle gaps 

The observed 2f energy of the neutron-rich Cr isotopes decreases monotonically 
with increasing the neutron number from Ei = 1007 keV in 56 Cr down to 446 keV in 
62 Cr, 2 ) suggesting a new region of deformation. The energy ratio E4/E2 between the 
2^ and 4^~ excited states becomes large E4/ E2=2.65 in 62 Cr 7 ) while the ratio is small 
E4/E2 ~ 2.2 for lighter isotopes. 6 ) The proton inelastic scattering experiment 7 -* in- 
dicates that the 2f state accompanies a large quadrupole deformation (3 ~ 0.25. 
Shell model analyses 2 )' 29 ) suggest a possible deformation for N >38 as an extended 

model space including at least the neutron I09 and 2ds orbits in the higher shell is 

2 2 

necessary to describe the decreasing trend of E2 with increasing the neutron number. 
There are several self-consistent mean-field calculations 21 ) -23 ) for the systematics of 
the ground state deformation, including the neutron-rich Cr isotopes. The calcula- 
tion using the Skyrme parameter set SkM* 21 ) produces a large prolate deformation 
around 62 Cr, but the other Skyrme models using SLy4, 21 ) SkP, 21 ) BSk2 22 ) and the 
relativistic mean-field model using NL-SH 23 ) do not. Here we analyse the microscopic 
deformation mechanism with focuses on the parameter set SkM*. 

First, we discuss the deformation energy curve in 56 ~ 68 Cr. We adopt SkM* and 
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Fig. 4. The quadrupole deformation energy curves of the neutron-rich Cr isotopes calculated with 
SkM*. The energy is arbitrarily shifted for the sake of comparison. 



the volume type pairing force (77 = in (2.1)) with the pairing strength Vq = -200 
MeV fm -3 . The value of Vo is adjusted to reproduce the odd-even mass difference 
in 56 ~ 66 Cr, which have the values between 1.5 and 2.0 MeV. The ground state values 
of (3, A n and A p are shown in Table II. The calculated deformation energy curves 
of the neutron-rich Cr isotopes are plotted in Fig. 4. It is seen that the deformation 
energy curve becomes soft toward the prolate direction with increasing the neutron 
number from iV=32. In 60 Cr the deformation energy is almost flat between /3=-0. 1 
and 0.3. In 62 > 64 > 66 Cr (iV=38-42), the minimum with a large quadrupole deformation 
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appears around (3=0.25 - 0.3, but then the deformation of the minimum decreases 
for N >42. The calculated isotopic tendency is consistent with the observed trend 
of Ei from iV=32 to iV=38. It is noted, however, that the energy difference between 
the spherical configuration(/3=0) and the quadrupole deformed minimum is small. 
In the case of 62 Cr, the difference is about 290 keV. If we estimate the zero-point 
energy in terms of the experimental value of £"2=446 keV, the deformation energy 
and the zero-point energy have comparable magnitudes. Thus it is suggested, within 
the SkM* model, that even in 62_64 Cr exhibiting the largely deformed minima a well 
developed static deformation is not realized, and that a large amplitude quadrupole 
motion of a transitional character is expected. The transitional nature seems to be 
consistent with the observed E4/E2 ratio 7 ) 2.65 in 62 Cr which lies in between the 
vibrator and the rotor limits 2 and 3.33. 

Next we analyse the mechanism of the deformation. For this purpose, we plot 
the Hartree-Fock single-particle energies as a function of (5. The single-particle orbits 
are obtained by rediagonalizing the HF Hamiltonian associated with the constrained 
Skyrme-HFB solution. The result is shown in Fig.5(a) and (b). The neutron subshell 



(a) SkM* neutron 



) SkM* proton 
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(c) SLy4 neutron 
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Fig. 5. The single-particle energies obtained with SkM*(left panels (a)(b)) and SLy4(right panels 
(c)(d)) for 62 Cr. The top panels (a)(c) are for neutrons, and the bottom panels (b)(d) are 
for protons. The thick solid line is the Fermi energy. The horizontal axis is the quadrupole 
deformation parameter (3. 



gap iV=40 between the /s and g$ orbits in the spherical region and an energy gap 

2 2 

iV=38 in the prolate region around /3=0.25-0.5 are notable. As the magnitude of 
the deformed iV=38 gap is comparable to the spherical iV=40 subshell gap, it is 
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easy to deduce that the deformed iV=38 gap may drive the nucleus toward prolate 
deformation with (5 ~0.3. We also note that the neutron gg orbit plays an important 

2 

role in forming the deformed states. Especially the Q = \ and | orbits which stem 
from the deformation splitting of vg$ are relevant as they exhibit the steepest down- 
sloping for (5 > 0. If neutrons occupy these vgi orbits, they cause a strong driving 

2 

force toward prolate deformation. This situation can be realized for (3 > 0.20. We 
show in Table II the occupation number for all of the vg<z_ orbits. The occupation 

2 

number ~3.0-3.8 in 62 Cr and 64 Cr, where the largest ground state deformation is 
realized, is consistent with the interpretation that the vgi fl = \ and | orbitals are 
largely occupied. As the neutron number goes beyond iV=38, the occupation number 
of vgg still increases slightly (3.8 in 64 Cr with iV=40 and 4.4 in 66 Cr with N=42). 
However, the quadrupole deformation then turns to decrease at N=42 since the 
upsloping orbits [303 |] and [301 \] stemming from vfs and vp\ are occupied. This 

explains the reason for the largest quadrupole deformation in 62 > 64 Cr with iV=38-40. 
The important role of the occupation of vg$ is consistent with the conclusions of the 

shell model analyses. 2 )' 3 )' 6 )' 29 ) 

3.2. Neutron-rich Ti and Fe isotopes 

It is interesting to investigate neutron-rich Fe and Ti isotopes around 62 ' 64 Cr 
with iV=38, 40 with which the largest deformation is realized in the Cr isotopes. 
Figure 6 shows the deformation energy curve calculated with SkM* for the iV=38 
isotones 60 Ti, 62 Cr, 64 Fe and for the iV=40 isotones 62 Ti, 64 Cr, 66 Fe. It is seen that 
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Fig. 6. The quadrupole deformation energy curves of the neutron-rich iV=38 and 40 isotones around 
Z=24 obtained with SkM*. 

the Cr and Fe isotopes have deformed minima at prolate deformation while the Ti 
isotopes do not exhibit this feature. Since the importance of proton configurations 
is easily inferred, we look into the proton Nilsson diagram in Fig.5(b), which is for 
62 Cr but essentially the same for the nuclei under discussion. We also show the 
occupation number in the proton /V orbits in Table II. The occupation numbers 

2 

are ~2, 4, and 6 in Ti, Cr and Fe respectively as naturally expected from Fig.5(b). 
Since the fl = \ and Q = | orbits stemming from tt/z have steep slopes in the 
prolate direction, simultaneous occupation of these orbits in Cr gives a large driving 
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Table II. The occupation numbers of the neutron 519 orbit and the proton fz orbit, the average 
neutron and proton pairing gaps A n , A p and the quadrupole deformation parameter f3 associated 
with the HFB solutions obtained with SkM*. 





neutron 09 


proton fz 


A n [MeV] 


A p [MeV] 


P 


60 Ti 


1.28 


2.02 


2.042 


1.198 


0.0 


b2 Ti 


1.28 


2.02 


2.096 


1.157 


0.0 


b« Cr 


0.355 


3.98 


1.619 


1.427 


0.0 


bU Cr 


0.794 


3.98 


1.970 


1.388 


0.0 


M Cr 


2.98 


4.00 


1.908 


0.405 


0.257 


b4 Cr 


3.85 


4.02 


1.904 


0.212 


0.262 


bb Cr 


4.39 


4.00 


1.833 


0.280 


0.247 


b «Cr 


5.40 


4.00 


1.756 


0.456 


0.222 


b4 Fe 


2.58 


5.99 


1.999 


0.280 


0.212 


bb Fe 


3.76 


6.00 


1.951 


0.008 


0.231 



force toward the prolate deformation. We can expect that the deformation is a less 
favored in Fe than in Cr because protons in Fe occupy the f2 = | orbits which is 
slightly up-sloping. Indeed the value of j3 at the minima in 64 ' 66 Fe is little smaller 
than those in 62 ' 64 Cr(cf. Table II , Fig. 6). On the other hand the spherical minima 
in Ti isotopes indicate that occupation of the Q = \ orbit alone is not enough to 
cause deformation. 

The observed 2f energies 30 ) in 64 ' 66 Fe are lower than that in 62 Fe, and the 
observed ratio 31) E 4 /E 2 =2.3Q in 64 Fe is smaller than that(=2.65) in 62 Cr. 7 ) The 
deformation energy curves indicating smaller quadrupole deformation in 64 ' 66 Fe than 
in 62 > 64 Cr are consistent with these experimental trends. In addition, the calculation 
predicts much less collectivity in Ti than in Cr and Fe isotopes. This is also consistent 
with higher E 2 in 58 Ti than in 60 Cr. 7 ) 

3.3. Relation to the N ~ Z =38-40 deformed region 

The role of the iV=38 deformed shell gap reminds us of the mechanism of the 
prolate deformation in proton-rich Sr and Zr isotopes with N ~ Z=38-40 as the 
large deformation in the proton-rich N=Z region originates from the presence of 
a deformed shell gap at N, Z=38 and 40 around j3 ~ 0.4. 26 )' 32 )' 33 ) It is therefore 
interesting to investigate how the deformation of the neutron-rich Cr, Fe nuclei is 
related to the N ~ Z=38-40 cases. Figure 7 compares the deformation energy curve 
in 62 Cr with those of 68 Zn(Z=30) and 76 Sr(Z=38) in the iV=38 isotone chain up 
to N=Z=38. Note that the deformation energy curve in the iV=Z=38 nucleus 76 Sr 
has a well-developed deformed minimum around (5 ~0.4. 

Figure 8 plots the neutron Nilsson diagram in 76 Sr. (The proton Nilsson diagram 
is not shown here as it is very similar to the neutron's except an overall energy shift 
due to the Coulomb potential.) Comparing the neutron Nilsson diagrams in Fig.8 
and in Fig. 5 (a), we see that the presence of the deformed iV=38 gap is common 
to both cases. The iV=38 deformed gap at [3 ~0.4 is universal in the sense that it 
exist both in the proton-rich and the neutron-rich regions. This confirms that the 
deformation mechanism in the neutron-rich Cr and Fe isotopes is intimately related 
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Fig. 7. The quadrupole deformation energy curves in the N=38 isotones obtained with SkM*. 
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Fig. 8. The neutron Nilsson diagram of 76 Sr obatained with SkM*. 



to that in the N=Z ~38 region. 

In spite of the above similarity, however, there is small but clear differences 
between the neutron Nilsson diagrams in the proton-rich N=Z ~38 region and in 
the neutron-rich Cr-Fe region. Looking closely at the N=38 deformed shell gap, 
we see that the relevant level crossing of the down-sloping vgv Q = | orbit and 
the up-sloping orbit from vfs and vpi occurs at smaller deformation (3 ~0.20 in 

2 2 

the neutron-rich 62 Cr than in 76 Sr, where the crossing takes place at [3 ~0.26. It is 
also seen that the iV=40 spherical gap at (3 ~0 is smaller in the neutron-rich 62 Cr 
than in 76 Sr. Both features are helpful to produce a largely deformed minimum 
at j3 ~0.27 in the deformation energy curve of 62 Cr even though the deformation 
driving effect due to the proton configuration is weaker in 62 Cr than in 76 Sr. Note 
also that the deformed shell gap at iV=40 around ~ 0.45 present in 76 Sr is barely 
seen in 62 Cr. Another difference is also seen in the iV=34 oblate gap which is 
large in the proton-rich region (cf. Fig. 8), leading to an oblate ground state in 68 Se 
(/3 = — 0.21). 13 )' 32 )' 34 ) However the iV=34 oblate gap is much weaker in the neutron- 



12 



H.Oba , M.Matsuo 



rich Cr isotopes(Fig.5(a)), and hence the deformation energy curve in 58 Cr(Fig.4) 
does not show very strong tendency toward the oblate deformation. 

3.4. Interaction dependence 

As we have seen above, the Skyrme-HFB model using the parameter set SkM* 
produces the onset of large quadrupole deformation in the N ~38-40 region of the 
Cr isotopes, being in qualitative agreement with the experimental observations. It 
is noted however that the calculations using other parameter sets such as SLy4, 
SkP 21 ) and BSk2 22 ) do not show largely deformed minima in the same isotopes. Let 
us discuss how the onset of deformation depends on the Skyrme parameter sets by 
comparing results which we obtain using SLy4 as well as SkM*. The calculation for 
SLy4 is the same as above except for the use of different Skyrme parameters. We 
show in Fig.9 the deformation energy curve in 62 Cr obtained with the parameter 
set SLy4. The deformation energy curve for SLy4 does not have a well deformed 




-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 
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Fig. 9. The deformation energy curve obtained with SLy4 in 62 Cr, plotted with the solid curve. 
The one for SkM* is also shown with the dotted curve for comparison. 



minimum, and it is stiff with respect to the quadrupole deformation. We show the 
neutron Nilsson diagram for SLy4 in Fig. 5 (c). By comparing with the Nilsson 
diagram for SkM* (Fig.5 (a)), we see a large difference in the position of the vgg 
orbit. The N=40 subshell gap at the spherical point is larger in the case of SLy4 
than SkM* by about 2 MeV. The higher position of ugg also shifts the position of 

2 

the 7V=38 deformed single-particle gap at larger deformation (3 > 0.35 in the SLy4 
case. These features apparently favor the stability of the spherical shape in contrast 
to the case of SkM*. All these observations indicate that the position of the gg orbit 
plays a central role for the onset of the quadrupole deformation in the neutron-rich 
Cr and Fe isotopes around iV=38. The precise position of the vgv orbit and the 

2 

size of the iV=40 subshell gap depend on the Skyrme parameter set. In other words, 
the deformation properties in the neutron-rich Cr region provides us with a rather 
strong constraint for a proper choice of the Skyrme parameter set. 
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§4. Conclusions 

We have developed an axially symmetric Skyrme-HFB code based on the 2D 
mesh representation in the cylindrical coordinate system in order to describe the 
quadrupole deformation of unstable nuclei. Our code has enough accuracy to analyse 
the deformation property as tested by comparisons with the calculations using the 
THO and the B-Spline methods. 

We applied the code to neutron-rich Cr, Fe and Ti isotopes around a possible 
new region of deformation with N >38. The quadrupole deformation energy curve 
obtained with the Skyrme parameter set SkM* indicates an onset of deformation 
in the isotopes with N ~38-42. The deformation energy curve is soft although the 
deformation of the minimum reaches (3 ~ 0.25 for iV=38 and 40, and hence the 
nuclei in this region are expected to exhibit strong quadrupole collectivity having 
transitional character rather than that associated with a well developed stable defor- 
mation. These results are in qualitative agreement with the trends of the presently 
available experimental data. The SkM* model describes also the Fe isotopes in the 
same mass region as transitional nuclei while the Ti isotopes are described as spher- 
ical nuclei having stiffer quadrupole deformation energy curve. By inspecting the 
neutron Nilsson diagram we discussed that the deformed iV=38 gap emerging at 
(3 > 0.25 as well as the position of vgg orbit play important role for the onset of 
quadrupole deformation. We have also shown that, contrary to SkM*, the param- 
eter set SLy4 is not able to reproduce the onset of deformation. The sensitivity to 
the Skyrme parameter set arises from the fact that the deformation in this region 
emerges as a consequence of a delicate competition between spherical and deformed 
configurations, for which the position of isgs orbit plays an essential role. 

2 
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